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Abstract: 

We compute the one-point probability distribution function of small-angle cosmic mi- 
crowave background temperature fluctuations due to curved cosmic (super-)strings with 
a simple model of string network by performing Monte Carlo simulations. Taking into 
account of the correlation between the curvature and the velocity of string segments, there 
appear non-Gaussian features, specifically non-Gaussian tails and a skewness, in the one- 
point pdf. The obtained sample skewness for the conventional field-theoretic cosmic strings 
is gi ~ —0.14, which is consistent with the result reported by Fraisse et al. Q. We also 
discuss the dependence of the pdf on the intercommuting probability. We find that the 
standard deviation of the Gaussian part increases and non-Gaussian features are suppressed 
as the intercommuting probability decreases. For sufficiently small intercommuting prob- 
ability, the skewness is given by < (a few) x 10 -2 . 
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1. Introduction 

Cosmic strings are line-like topological defects formed in the early universe through a 
spontaneous symmetry breakdown. It was claimed that the formation of cosmic strings 
at the end of inflation is a generic feature of supersymmetric grand unified theories |2j. 
The string tension [i is directly related to the symmetry breaking energy scale. Also 
recent studies of stringy cosmology have revived interest in cosmic strings, because it was 
pointed out that brane inflation models may produce another class of string objects, called 
cosmic superstrings 0, |4|, [5|, [6[. Cosmic superstrings may have different properties from 
those of conventional field-theoretic cosmic strings. One of the observationally interesting 
differences is the intercommuting probability P. It can be significantly smaller than unity 
for cosmic superstrings , while normally P = 1 for field-theoretic strings || (but see Q ) . 

The imprint of cosmic strings on the cosmic microwave background (CMB) has been 
widely studied. Although cosmic strings were excluded as a dominant source of the observed 
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large angular scale anisotropy [10, 11, [T^] (see also |H| for a review), a signal due to 



cosmic strings could still be observed at small angular scales jj], [l(J with future arcminutes 
experiments such as South Pole Telescope [14] or Atacama Cosmology Telescope [15]. On 
small angular scales, the primary fluctuations are damped and only the integrated Sachs- 
Wolfe (ISW) effect is relevant. The Gott-Kaiser-Stebbins (GKS) effect |§, is the most 



characteristic signal of the ISW effect due to cosmic strings. The GKS effect is caused 
by a discontinuity of the gravitational potential across a moving string segment when 
photons passes by the moving string. If photons are scattered by a number of moving 
string segments, the observed temperature fluctuations appear as a superposition of the 
discontinuities. 

Since the network of strings is a highly nonlinear object, its non-Gaussian features 
may help us distinguish cosmic string signals from other secondary effects and hence may 
enhance the observability. A skewness, which is a measure of the asymmetry of a probability 
distribution function, is one of the simplest features of non-Gaussianity [18]. The skewness 
of temperature fluctuations is defined by 



where A = (T — T)/T = AT/T, the bar " ~ " denotes the statistical average over a CMB 
map and cta the standard deviation. Recently, Fraisse et al. [0 found that the one-point 
probability distribution function (pdf) of the temperature fluctuations due to conventional 
cosmic strings has non-Gaussian tails and a negative skewness, g\ ~ —0.23. 

Non-Gaussian features would also appear in the bispectrum [19]. Some authors [20, 



21, |22|] calculated the bispectrum and the trispectrum by using the formula derived by 



Hindmarsh [23|. In particular, Hindmarsh et al. |20| pointed out that the numerically 
measured amplitude corresponds to a |/nl| ~ 10 3 for the local-type of primordial non- 
Gaussianity. 

In pi| , we computed analytically the one-point pdf of small-scale temperature fluctua- 
tions with a simple model of long straight segments and kinks. We derived basic equations 
for the evolution of string segments in cosmological backgrounds by extending the velocity- 



dependent one-scale model with the intercommuting probability P / 1 [|24| , [25 , 26| , [27]] . It 
was found that the obtained one-point pdf consists of a Gaussian component due to fre- 
quent scatterings by long straight segments and non-Gaussian tails due to close encounters 
with kinks. The dispersion of the Gaussian component is consistent with the result in Q 
and the non-Gaussian tails also can be fitted to that in [Q] by using two phenomenological 
parameters. The effect of the intercommuting probability P was also studied there. It 
was shown that the non-Gaussian tails diminish as P decreases. However the obtained 
one-point pdf is symmetric for positive and negative temperature fluctuations, hence did 
not reproduce the non-zero skewness reported in jlj. 

In this paper, we study in depth the non-Gaussian features, especially skewness, of 
the CMB temperature fluctuations due to cosmic (super-)strings, and interpret numerical 
results obtained in [|J. For this purpose, we adopt a simple model of string network and 
perform Monte Carlo simulations. We consider the correlation between the curvature and 
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the velocity of string segments as the origin of skewness and study the dependence of the 
skewness on the inter commuting probability P. 

This paper is organized as follows. In section 0, we give a brief review of the derivation 
of basic equations for the evolution of string segments incorporating the intercommuting 
probability P [24|. In section |3| we introduce general formulae for the temperature fluc- 



tuations due to cosmic strings ^8|, |2!j. We also discuss the correlation between the 
curvature and velocity of string segments which gives rise to the skewness in the one point 
pdf. Then, in section ^j, we give the result of our Monte Carlo simulations. We show that 
the obtained one-point pdf has a negative skewness when the correlations are taken into 
account. The skewness is plotted as a function of P. Finally, we summarize our results in 
section [|. 

2. An analytic model of the string network 

A string worldsheet can be described by = x^(a a ), where x^ and a a are the spacetime 
coordinates and worldsheet coordinates, respectively. The induced metric of the worldsheet 
lab is given by 

dx^ dx v 

la b = 9,^^, (2.1) 

where is the spacetime metric. We assume that the action for the string dynamics is 
well approximated by the Nambu-Goto action: 

j2 



Sng = M / V=ld 2 a , (2.2) 



where [i is the tension of the cosmic string. 

Let us consider string dynamics in Friedmann-Lemaitre-Robertson- Walker universes 
with the metric 

<*4lhw = « 2 (^?) ("<V + dr2 ) ■ ( 2 - 3 ) 
We choose the temporal gauge : 

a = v , a 1 =a, r V = 0, (2.4) 

where the bold letters denote the 3-vectors on the comoving space and the dot and the 
prime denote the derivative with respect to rj and a, respectively. The equations of motion 
on the cosmological backgrounds are given by |30| ] 

f + 2%{l-r 2 )r = \{^\ > (2-5) 
e = -2Her 2 , (2.6) 
where % = a/a and e is the energy per unit coordinate length defined by 



r /2 "1 1/2 



(2.7) 
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We define the total energy E and the average root-mean-square velocity of a string: 

E = pair,) J eda , v 2 ms = (r 2 ) = • (2-8) 

In the velocity-dependent one-scale model (VOS), the string network is characterized by 
only two physical quantities: the correlation length £ and the rms velocity t? rms . The 
correlation length £ is defined by 

V p 

where p is the total string energy density, which related to E as 

pocE/a 3 . (2.10) 

In our treatment, we also take the energy loss due to loop formations into account. A 
loop formation can occur through the intercommutation of two segments or self-intercommutation 
of a single segment. The characteristic time scale of the interval of loop formations with 
the intercommutation probability P is ~ £/(Pv Tms ). Then the energy loss due to the loop 
formation can be described as 

= -cPv? , (2.11) 

at / to loops ? 

where we have introduced c as a constant which represents the efficiency of loop formations. 

For brevity, we use the quantity 7 = 1/H£ rather than £ with H = aH. For a universe 
with the scale factor a(t) oc t" with physical time t = J a(j])dr), the equations of motion 
for 7 and i> rm s can be derived from Eqs. Q - (pd|) as [||, ||, |f| || 

- d J7 = ^-P-\PcPVr m sl-Pv 2 ms , (2.12) 

7 at 2 



™rms / -1 9 
- (1 - V* t 



k 

2Hv r 

R 



(2.13) 



dt imsJ 

where we have replaced ((^ 2 ) 2 ) by (r 2 ) 2 . The momentum parameter k and the curvature 
radius R are defined below. For the curvature term, we set 

where we have introduced the proper coordinate length measure along a string by ds 
t' (in = vT — r 2 eda and u is the unit vector parallel to the curvature radius vector on the 
comoving space. The momentum parameter k is defined by 

k ^ ((l-r 2 )(r.u)) = R L.l(lL)'\ ( 2 15 ) 

In this paper, we use an approximated form of k |^5| : 

M*W) = ^^K (2.16) 
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It is known that a string network approaches the scaling regime where the characteristic 



scale grows with the horizon size [31, p^j. For our Monte Carlo simulations, we assume 
that the scaling is already realized by the time photons leave the last scattering surface. 
This means that 7 and t> rms are constant in time. Hence we set the left-hand sides of 
Eqs. ( gig ) and ( gig ) equal to zero. We solve them numerically with the assumptions that 
the universe is in the matter-dominated era, /3 = 2/3, and the curvature radius is equal to 
the correlation length, R = £. 



3. CMB Temperature Fluctuations 

3.1 The Hindmarsh-Stebbins-Veeraraghavan formula 

A formula for the temperature fluctuation due to a Nambu-Goto string was derived by 



Stebbins and Veeraraghavan [28, ^] and by Hindmarsh 22]. The Hindmarsh-Stebbins- 



Veeraraghavan (HSV) formula is given by 

AT f (1 + n • f) [(In + X) ■ u] 



AT f 1 

A(n , r obs , 7? obs ) = — = -2G/x J da- 



(X + n ■ X) (X - r ■ X) 



, (3.1) 



where X(a) = r Q ^ s — r(a,rji c (cr)) is the comoving position of the observer relative to that 
of the string, X = \X\, u is defined by 

/ n ■ r' \ , . 

u = r - r r' , 3.2 

\l + n-r) v ; 

and ^ic(c) is the conformal time along the intersection of the observer's past light-cone and 
the string worldsheet, 

%bs - ?7ic(o") = X(a) . (3.3) 

It should be stressed that the temperature fluctuation depends on the string distribution 
on the observer's past light-cone. 

The integration range £ must be appropriately determined. In our calculation, we 
consider only a segment of a long cosmic string at each scattering. Namely we take £ as 

E : X ± (a) < (3.4) 
a 

where 

X ± (a) = X(a)-(X(a)-n)n. (3.5) 



This procedure is consistent with the treatment in our previous paper [24|. 

It is useful to adopt the small angle approximation to understand various effects con- 
tained in the HSV formula. In the small-angle approximation, we assume that the angle 
between the direction of the source —X and the line-of-sight n is sufficiently small. This 
situation is realized when the length of a string segment, which is approximately equal to 



- 5 - 



the Hubble radius at the epoch of scattering, is much smaller than the distance between 
the string and the observer. Thus we have 

\ 2 

(3.6) 



A <1, X(a) sa -In + X ± + n + 



X(a) 



2X 



This approximation is valid as far as a scattering due to a segment at large redshift z>l 
is considered. 

Let us consider the leading order in the small angle approximation, 



-4G[i / da 



X 1 - ■ u 



(X^f 



V=Vlc((7) 



For convenience, we decompose this into two parts as [28] 
A « A" + A x , 



A 11 = AGfj, / da a 



X 



{X^f 



A ± = AG/jl I da 



a± 



mc(o-) 



where the coefficients an and a± are defined by 



(3.7) 

(3.8) 
(3.9) 

(3.10) 



mc(o-) 



o II 



n ■ r r r 

— — ^ , a± = n • — 7 x — : 



(3.11) 



Note that d/dcr in the above formulae is the derivative along the intersection of the past 
light-cone with the worldsheet, 



d_ _ rf7^ c (o-) 3 



We may rewrite Eq. ( |3.9| ) by integration by part as 

x 1 - 



dan 

-AGfi / dcr— —In 
da 



L 



cut 



(3.12) 



(3.13) 



^lc(o-) 



where L cut is the cutoff scale of the scattering, L cut = £/a. It is easy to see from this 
expression that for an exactly straight and uniformly moving segment, that is for r = const, 
and r' = const., A" vanishes. 

If the segment has a kink at a = a\- , we have 

aj| = aj +) e(a - cr fc ) + aj _) G(a & - <r) , (3.14) 

where aj| are the constant parameters characterizing the kink and 0(c) denotes the 
Heaviside step function. Then, reproduces the kink temperature fluctuation, 

X ± (a k ) 



A 11 « A kink = -4G/ia kink In • 



-C'kink 



(3.15) 
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where we have defined akink = — oci. ' and Lkink is a mean distance between kinks 
1 24 , ^9], 28]. If the curvature of the segment is taken into account, dan/da 7^ and All 
appears. Therefore, A" can be regarded as contributions of the deficits on the segment 
(kinks and cusps) and the curvature of the segment. 

For an exactly straight and uniformly moving segment, the position of the segments, 
X- L (a), can be written as 



.,(-) 



X ± (a) = S + 



a 



dX- 



da 



(3.16) 



where we have introduced the impact parameter 8 = X^(a = 0), 5 = \8\, and the unit 
tangent vector e oc dX 1 - /da. Then A 1 - becomes 



A 



-4Gfi 



VT 



arctan 



dx 1 - 



\a + 8 cos tp 



8 sin <p 



VT 



-a±_ arctan 



\ad sin tp I 



(3.17) 



where 8 ■ e = 5 cos (p. It is clear to see that the dominant contributions are given from 
the string portions with \a\ < 6/\ v. I. This is just the Gott-Kaiser-Stebbins (GKS) 



effect [16, 17]. Hence A 1 - may be regarded as the generalization of the GKS effect. If the 
curvature of the segment is taken into account, correction terms due to higher derivatives 
will appear in A -1 -. 

Although most scatterings occur at z > 1, some photons may be scattered by a 
segment at z ~ 0{1). Therefore it is desirable to check the validity of the small angle 
approximation. For this purpose, in our Monte Carlo simulations, we have taken into 
account the next order corrections and compared the result with the one obtained in the 
leading order approximation. We found no difference between these two, supporting the 
validity of the small angle approximation. 

3.2 The Gott-Kaiser-Stebbins effect 



Next we consider the GKS effect [16, 17 1. Using Eq. ( 3.17 ), the rms temperature fluctuation 
due to the GKS effect of a string is estimated as 



^GKS 



a scg Gji arctan ( -^r 



v 



rms 



0.8 



(3.18) 



where ct s 



„2 \V2 



o" ) . In the case of a close encounter, i.e. 8 <C £/a, this reduces to 



A 



GKS 



47T- 



y/1- r 



2 

rms 



(3.19) 



Thus, this produces the discontinuity of the temperature fluctuations across the moving 
string segment when photons passes by the moving string. 

The temperature fluctuation due to a segment approaches zero when the impact pa- 
rameter is larger than the segment length £. This implies that a segment has an effective 
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Figure 1: The standard deviation of the Gaussian components of the pdf due to GKS effects, ctgks 
as a function of intercommuting probability P. 



cross section ~ £ 2 . In terms of £, the number of segments in a horizon volume is expressed 
as N seg = l/^H 3 = 7 3 . Therefore the optical depth for a CMB photon is 



'opt 



~LHS 



3 e 2 



dz 



16 



seg * H(l + z) 
cP N - 1/2 



7log(l + zlss) 



0.23 



(3.20) 



where 7 « (3C-P/V2-7T 2 ) -1 / 2 , and we have put zlss ~ 1100 and used c « 0.23 as the 
standard value [j33"|| . Since r opt 3> 1, a photon ray is scattered by segments many times 
on its way from the last scattering surface to the observer. Therefore, the temperature 
fluctuations behave like a random walk and the one-point pdf becomes well approximated 
by a Gaussian distribution as a result of the central limit theorem. Thus the GKS effect 
gives the Gaussian component of the one-point pdf, 

dP 

dA GKS 



1 

/ 27TCTGKS 

where the standard deviation is estimated as 



-A 2 /24 KS 



(3.21) 



CGKS 



A 



GKS" 



13^ 



/ Q seg \ 



CP 



-1/4 



(3.22) 



Vl/^27 VO-23^ 

which is shown in Fig. |l| as a function of intercommuting probability P. In [24], we found 
that the above estimate of the standard deviation agrees well with the result of numerical 
simulations (for P = 1) by Fraisse et al. 0. 

3.3 Correlations as sources of the skewness 

The one-point pdf of the temperature fluctuations due to the GKS effect is symmetric under 



A —7- — A since the form of the GKS effect is obviously symmetric as is seen from Eq. (3.1 
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However, a skewness may appear if we take account of the curvature of segments and its 



correlations with the velocity, as discussed in [24, 20]. If a segment has a curvature, the 
probability of producing positive or negative temperature deviation becomes asymmetric. 
The amplitude of the asymmetry depends on the angle between the curvature vector and 
the velocity vector of the segment. However, if the curvature and the velocity are not 
correlated, the statistical averaging results in a vanishing skewness. Therefore, a correlation 
between the velocity vector and the curvature vector is necessary to produce a skewness. 
In this paper, we investigate two kinds of correlation. One is due to the breaking of the 
time reversal symmetry caused by the cosmic expansion, which has been already pointed 



out in [20[. The other is due to the effective curvature of a string on its intersection with 
the observer's past light-cone. We call this "light-cone effect". 



3.3.1 Correlation due to cosmic expansion 

From Eq. (§T| ), assuming the scaling regime, we have 



A: 
R 



av r 



1 



2Hv r 



(3.23) 



It is easy to see that a nontrivial correlation between the velocity vector r and the curvature 
vector r" exists if there is Hubble expansion. Each configuration of string segments must be 
consistent with this correlation. We use this expression when we perform our Monte Carlo 
simulation. A method for determining string segment configurations in our simulations is 
explained in section 13. 



3.3.2 Light-cone effect 



In the HSV formula, Eq. (3.1), the integration is along the light-cone. In the small angle 
approximation, the coordinates of a string on the past light-cone emanating from the 
observer, X(a), can be expanded as 

v , \ v (dX\ 1 (d 2 X\ 2 ml , 
X(a) = X +( — J a+-[ — \ n- + 0(a> 



Xo 



da 2 



r'o + 



di] 



da 



da 

2 



ro 



a 



ro + 2 



o 



da 



d 2 y\c 
da 2 



a 2 + 0(a A 



where dr]\ c /da and d rj\ c /da are estimated from Eqs. ( |3.3D and (]3.6|) as 



drjic 
da 

d 2 y\c 
da 2 




(3.24) 



(3.25) 



+n 



(3.26) 
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It is remarkable that even if a string is exactly straight and moves uniformly, r$ 
tq = Tq = 0, it always has an "effective curvature" . Namely, 

ffr\ _ ffrn. 



where note that d/da is a derivative along the light-cone; see Eq. ( |3,12| ), This obviously 
means the existence of a correlation between the effective curvature and the velocity of the 
segment. However, this effect is subdominant in the small angle approximation, because of 
the suppression factor a / Xq in Eq. ( [3.24 ) , where Xq is of the order of the angular diameter 



distance from the string to the observer. Thus as far as scatterings at z > 1 are concerned, 
the light-cone effect is smaller than the effect due to cosmic expansion. 

4. Monte Carlo simulations for the one-point pdf 
4.1 Setup 

We compute the one-point pdf of the small-angle CMB temperature fluctuations due to 
curved string segments with a simple model of the string network (VOS) by performing 
Monte Carlo simulations. 

4.1.1 Scattering probability 

Curved segments are assumed to be located randomly between the LSS and the present 
time consistently with the VOS model. We divide the redshift range < z < zlss into a 
number of bins with the width 

A log(l + z) = lQg(1 + ZLSs) Ap = , (4.1) 

where Ap is the scattering probability assigned to each bin. The scattering events due to 
cosmic strings can be randomly arranged on the observer's line-of-sight by using Eq. ([4.1|). 
We set Ap = 10~ 2 in our simulations. 

4.1.2 Assumptions and the expansion of the HSV formula 

To take into account the curvature of a string segment, we expand X- L (a) to second order 
m a as given by (p|). We assume the higher order terms are negligible. Hence we set 



^-' t ®I + '(^' (4.2) 



X X V da J X 2 V da 2 J X 

We also assume the validity of the small angle approximation and assume S/Xq = 0(e) 
and \ d ^ a \q(j/Xq = 0(e), where e is a small expansion parameter. Furthermore, we assume 
that the curvature r'^a , the rotation r'^a and the acceleration r^a of a string segment can 
be treated as perturbations from a uniformly moving straight string segment, 

\r%\cr ~ \r \a ~ \r' \cr = O(e) <^l. (4.3) 
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The above order estimate comes from the fact that r"a ~ er/£ ~ ct/Xq. Note that these 
order estimates imply that the first two, leading order terms in (|4.2|) are already of O(e) 
and hence the last curvature term, which is 0(e 2 ), is only 0(e) relative to the leading order 
terms. 

Now we can expand the reduced HSV formula Eq. ( |3.7[ ). As discussed in the above, 
the effect of the curvature appears 0(e) relative to the leading order terms. Therefore, we 
expand the numerator X 1 - ■ u to 0(a 2 ) and the denominator (X^) 2 to 0(cr 3 ). Hence it 
can be written as 



Ei=o a i ai 



A « -AGfi I da ^=° n l \ , (4.4) 



where a\ (i = 0, 1, 2) and j3j (j = 0, 1, 2, 3) are given in Appendix |A|. Once the string con- 
figuration is fixed at each scattering event, we can immediately calculate the temperature 



deviation using Eq.(4.4). 



4.1.3 Degrees of freedom of a string configuration 

Our remaining task is to determine string configurations. A string configuration can be 
characterized by several random variables. Details about the random variables are discussed 
in Appendix |B[ Here, we just count the number of degrees of freedom (dofs) to determine 
a string configuration. 

For a string segment with the curvature taken into account, originally we have 3 dof for 
t*o, 6 dof for fo and r' , and 9 dof for fo, r' and r '. For ro, using the coordinate dof, there 
remains only 1 dof, which is the impact parameter 5. For fo and r' , the number of dof is 
reduced to 4 by using 2 constraint equations (the gauge condition and the normalization) 
for the string dynamics. The remaining dof are the magnitude of the velocity v and the 3 
angular parameters {0, (ft, ift) that determine the direction along the string and the velocity 
vector in the plane orthogonal to the string direction. Finally, the number of dof of fo, 
r' and r ' is reduced to 4 using 3 equations of motion and 2 constraints. One of them 
corresponds to the inner product of the velocity vector and the curvature vector, fo • r '. 
The remaing 3 dof correspond to the curvature vector orthogonal to f o and the rotation 
of the string segment orthogonal to both fo and r '. Lacking the knowledge of how much 
freely moving string segments would be curved, we set these remaing 3 dof to zero for 
simplicity. Thus we have 6 dof in total, (5, v, 6, (ft, ip, fo • t'q). 

In this paper, we further reduce the dof by fixing the velocity and using the scaling 
assumption. First we set v = v rrns . In the scaling regime, the curvature of a string segment 
and the velocity must be correlated. Locally at each event of scattering, one may neglect 
the Hubble expansion, and we have 

1 / 1 / r' \ \ / f • r"\ fn-r'n 9 , > 

' - ' < \ - / ^ - 11 ~ 2tf< s (1 - v 2 ms ) . (4.5) 



-< r 

a \ e 



Then 

f • r'o = 2aHv 2 ms (l - v 2 ms ) . (4.6) 

Thus, the remaining dof are (5, 8, (ft, tft). For a demonstration of the temperature fluctuation 
due to a curved segment, we show the CMB map with 6 = (ft = tt/2 and ift = in Fig. || 
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Figure 2: The amplitude of the curvature vector normalized by £ = I/-H7, |j"o'|£ = (r ' •r'o/ u rms)£ 
0(e) as a function of P. 
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Figure 3: The map of the temperature fluctuations due to a curved string segment with 9 = <f> = 
7r/2 and tp = 0. Also we define (Hx,Hy) as the Hubble normalized two-dimensional coordinate 
relative to cr = in the small patch of sky. The green line denotes the position of the string segment, 
X (cr), the dashed lines are contour lines with the width lG/i and the red vectors are the velocity 
vectors at each point. 



4.2 Numerical Results 

We performed Monte Carlo simulations with 10 7 photons on NEC Express 5800 XeonMP 
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Figure 4: The one-point pdf of the CMB temperature fluctuations induced by curved cosmic 
strings and kinks for the intercommuting probability P = 1. The solid curve (in red) is the one- 
point pdf obtained by averaging over scattering of 10 7 photons by curved string segments, added 
on top of an analytic pdf for kinks obtained in |24|] . The dashed line (in blue) is the best Gaussian 
fit. The dotted line (in green) represents the one-point pdf for straight segments and kinks. 



3.16GHz at Yukawa Institute for Theoretical Physics, Kyoto University. It took two CPU 
days for P = 1, three CPU days for P = 0.5 and five CPU days for P = 0.1. 

For P = 1, the obtained one-point pdf of the total temperature fluctuations is shown 
in Fig. In order to take into account the contribution from kinks, we just added an 
analytic pdf for kinks obtained in our previous paper (|1|). The one-point pdf significantly 



deviates from the best Gaussian fit. The sample standard deviation is cx s i m 

(p = 1) « UG/x, 

consistent with the analytical estimate in [^4|] and Fig. |l| as well as with the numerical result 
of M. The sample skewness measured in Fig. is 



(A - A) 3 

gi (P=l)= { - ^L«- .14. (4.7) 

Thus the one-point pdf has a negative skewness. 

In order to investigate the dependence of the one-point pdf on the intercommuting 
probability P, we computed the one-point pdfs for P = 0.5 and P = 0.1 and the obtained 
one-point pdfs are shown in Figs. [5] and |6|. The standard deviations are cr s ; m 

(P = 0.5) » 

18G/i for P = 0.5 and er s j m (P = 0.1) ~ 31Gfi for P = 0.1. These are consistent with the 
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Figure 5: The same as Fig. 4 but for the intercommuting probability P = 0.5. 



analytical estimate in [24| and Fig. |]. The sample skewness was found as 

9l{P<1) = {- 0nt °< P = *' (4.8) 
y V ' I -0.04 for P = 0.1. V 7 

Furthermore, we calculated the skewness for various values of P and obtained a power-series 
best fit in Fig. [7|. As we see, the dispersion of the Gaussian part increases as P decreases 
(see Fig. ||) and the non-Gaussian features, i.e., the skewness and the non-Gaussian tail, 
are suppressed. Our result suggests that the sample skewness for a sufficiently small P 
approaches g\ < (a few) x 10~ 2 . 



5. Summary 

In this paper, we have computed the one-point pdf of small-angle CMB temperature fluc- 
tuations due to the "curved" cosmic (super-)string segments by performing Monte Carlo 
simulations. Our purpose was to examine the numerical result obtained in Q with a 
simple model of the string network under a few reasonable assumptions, i.e., the velocity- 
dependent one-scale (VOS) model. We have found that the presence of the nonzero mo- 
mentum parameter k corresponds to the correlation between the curvature and the velocity 
of segments, and therefore the Hubble expansion is essential for a nontrivial correlation (see 
Eq. ( 2.15| ) and Eq. Q3.23 )). We have also found another type of correlation between the 



curvature and the velocity induced by the effective curvature along the intersection of a 
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Figure 6: The same as Fig. 4 but for the intercommuting probability P = 0.1. 



string worldsheet and the observer's past light-cone ("light-cone effect"). The light-cone 
effect always exists independent from the properties of string configuration and the cosmic 
volume expansion. 

Combining these two effects, we have calculated the one-point pdf numerically and 
showed that it reproduces features of the numerical simulations in ]lj very well. In partic- 
ular, we have found that the one-point pdf has a negative skewness. It was also found that 
as P decreases the standard deviation of the Gaussian part increases and the non-Gaussian 
tails and the skewness are suppressed. Our result suggests that, for the sufficiently small 
P, the skewness approaches < (a few) x 10 -2 . 

Note that the observed temperature fluctuations are given by contributions not only 
from cosmic (super-)strings but also from primary and other secondary anisotropies. It 
is known that conventional cosmic strings can contribute at most 10% of the total power 
at I = 10 in the CMB fluctuations [jllj. Thus, in real world, contributions from primary 
anisotropies make the skewness significantly small on large angular scales. In order to 
suppress this effect of the large-scale anisotropies, it would be effective to choose an obser- 
vation area of a sufficiently small angular size (~ O(10)arcmin 2 ) where the primary CMB 
fluctuations have been completely damped out. This would remove the contribution of the 
primary fluctuations to the denominator of the skewness formula (see Eq. ( |1.1[) ), hence 
enhance the observability of the skewness. 

It is known that O(0. 1) skewness would indicate nonlinear parameter |/nl| ~ O(10 3 ) f34| , 



2C]. In this paper, we have shown that, for smaller intercommuting probability P, non- 
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Figure 7: The skewness as a function of the intercommuiting probability P. The solid dots (in 
red) are the skewness obtained by simulations, each with 10 photons. The dashed line (in blue) is 
the power-series best fit function: gx(P) « 0.11P 2 - 0.23P - 0.02. 

Gaussian features become more difficult to find. To discuss the relation between the skew- 
ness and the non-Gaussian parameter /nl m the bispectrum, it is important to compute 
the angular power spectrum in our model. The work along this direction is in progress |J5|] . 
Since the non-Gaussian features significantly depend on the intercommuting probability, 
they could help us to distinguish field theoretic cosmic strings and cosmic superstrings. 
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Figure 8: A schematic picture of the surface perpendicular to the line of sight n at scattering. 




A. Reduction of the HSV formula 



The reduced HSV formula is 



A = -4Gfi / da- 



u 



1 



X 



(X 1 ? 



(A.l) 



V=V\c(cr) 



where we have introduced u = (1 + n ■ r)u, and added terms second order in the small 
angle approximation. 
We set 



X ± (a) = 6 + 



dX 1 
da 



(A.2) 



* \ uu / 

where we have introduced the expansion parameter e such that 5/Xq = 0(e), etc., as 



discussed in 4.1.2. The absolute value of X is then expanded as 



5 2 + 2 



S ■ 



dX J 
da 



a + 



+ 



dX- 



o 

d 2 X ± 



dx- 



da J V, da 2 j „ 



da 

a' + O^X 2 



2 {S U 2 x- 

o V da 2 /0 



(7 



= a + aia + a 2 cr 2 + a 3 cr 3 + 0(e 4 A^) . 

We also have another useful expression, 

X X^ • r 
I- — .r^l + nr — = L + L 1( r + 0{e 2 ) , 



(A.3) 



(A.4) 
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where 



L = 1 + n • r 

1 , ! n-— ^ — 

"' daj X 



o ■ r 
1 



da 



• T 



Then the HSV formula due to a curved segment can be approximated as 



A = -4Gfi / da 



2^=0 a i° 



2 E^o/^' 



where the coefficients of the numerator are 



a = S ■ u + 



2X n 



(n • u ) , 



"2 



.<Wo 



dX 1 
do" 



«o + 



1 







V d(T /o \ da )o 2 



2X 



ai(n • tto) + ao 



" 



• w 



and those of the denominator are 



/3 = L ao , 
/3i = L ai + Lia , 
/3 2 = L a 2 + iiai , 
fa = L a 3 + Lia 2 . 



(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 
(A.10) 



(A.ll) 
(A.12) 
(A.13) 
(A.14) 



A schematic picture of the relation between the normalized bases and the origin of the 
worldsheet coordinates is given in Fig. |8| 



B. Random Variables 

For our Monte Carlo simulations, we need to specify random variables in a physically 
reasonable way. In our situation, the random variables are those which specify the string 
configuration at each scattering event. They are the string position tq, the velocity vq, the 
direction r' , the curvature r ', the rotation r and the acceleration r '. Thus there are 18 
parameters in total. However, not all of them are independent from each other because of 
the equations of motion and constraint equations. 

Since at each scattering event, the cosmic expansion can be neglected, the equations 
of motion can be approximated by those on the Minkowski background. In the conformal 
and temporal gauge, the equations of motion are given by 

r-r" = 0. (B.l) 

The gauge conditions lead constraint equations as 

f V = 0, r 2 + r' 2 = 1. (B.2) 
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Differentiating the constraint equations, we have 



r'o ■ r = -r' Q -r' , r • r' = -r ■ r' . 



Using Eqs. ( |B.1| )-( |bT3|) , we can derive the following expressions: 



r' = yl - v 2 n» , t-q = vt\ 



VI — v- 



= K-*l)»[|+W-*l)*l + W-*2)*2, 



vT 



(B.3) 

(B.4) 
(B.5) 

(B.6) 



where (n\\,ti,t2) are unit basis vectors with nu parallel to r' and t\ parallel to tq. In 
terms of coordinate basis vectors (n, si, S2) with n being the unit vector along the light of 
sight, those parameters can be written as 



v = \r \ , 

nil = sin cos (ps\ + sin sin (pS2 + cos On , 
ti = (— sin ip cos cos <f> — cos ip sin 0)si 

+(— sin ip cos sin <p + cos ^ cos (j))s2 + sin ^ sin On , 
<2 = (— cos ip cos cos <p + sin ^ sin 4>)sx 

+(— cos V' cos sin — sin ip cos </>) s 2 + cos ip sin #n , 

where the ranges of the angular coordinates are 

< <9 < vr , < </> < 2vr , 0<V<2vr. 



(B.7) 
(B.8) 

(B.9) 

(B.10) 

(B.ll) 



Thus the total number of degrees of freedom (dof) is reduced to 9; (5, v, (f' ■ ti), (Tq • ti), 
(r' -t 2 ),(r%-t 2 ) : e,cP,iP). 

The impact parameter 6 = Xq has to be chosen randomly at each scattering. The 
impact parameter has two dof. One is the absolute value 6 and another is the angular 
coordinate Ox- Then 



X =5 cos OxS\ + sin6*xS2 
Performing a 2-dimensional coordinate transformation, 



si 

S2 



cosx smx 
— sin x cos x 



81 



(B.12) 



(B.13) 



the impact parameter Xq can by expressed as 



X^ = 5 



cos(O x + x)h + sin(6*x + x)$2 



(B.14) 



Therefore, by using the dof for the rotation of the 2-dimensional coordinate basis the 
angular dof Ox can be chosen freely. We set Ox = vr/2. 
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For the distribution of 5, from its definition we set P(5)d5 oc 5d5 with the cutoff at 
5 = £/a. Then we fix the range of integration X ± (a) < £/a for a given 5. For a string 
segment with the parabola-type shape, Eq. (fL2]), the condition can be written as 

ao 7T + aicr + 02<t 2 + a3<r 3 < . (B.15) 

gt 

For a given value of the impact parameter at 5 < £/a, the range of a is chosen to satisfy 
the above equation. 

For the angular variables 9, <f> and ip, we choose uniform distributions P(cos 9) = 
1/2 , P(4>) = P(ip) = 1/2-7T. As for v, we fix it at v = v rras for simplicity. 

We fix the value of {t'q -t\) by using Eq. (|4.5| ) so that it is consistent with the correlation 
between the string velocity and the curvature. As for the remaining three dof, we set 
( r o ■ t%) = 0, (r' ■ t\) = and (r' ■ (2) = 0. Note that this is solely due to the lack of our 
knowledge about the curvature and rotation of string segments. 
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